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SUMMARY 


Several  variational  problems  involving  optimum  wing  and  body  combina¬ 
tions  having  minimum  wave  drags  for  different  kinds  of  geometrical  restraints 
are  analyzed.  Particular  attention  is  paid  to  the  effect  on  the  wave  drag 
of  shortening  the  fuselage  and,  for  slender  axially  symmetric  bodies,  the 
effect  of  fixing  the  fuselage  diameter  at  several  points  or  even  of  fixing 
whole  portions  of  its  shape . 


INTRODUCTION 


Recently  several  authors  have  used  linearized  theory  to  study  the 
wave  drag  of  wing -body  combinations  traveling  at  supersonic  speeds  (see, 
e.g.,  refs.  1  to  5) •  These  studies  have  clearly  demonstrated  the  impor¬ 
tance  of  finding  the  wave  drag  of  a  whole  airplane  rather  than  the  separate 
wave  drags  of  its  various  parts  (wings,  fuselages,  etc.),  since  the  magni¬ 
tude  of  the  interference  terms  can  predominate.  In  effect,  this  means 
that  various  optimization  problems  for  bodies  -  such  as  the  problem  of 
finding  the  body  shape  having  a  minimum  wave  drag  for  a  given  volume  - 
should  be  re-examined  when  interfering  wings  or  other  bodies  are  in  the 
same  flow  field.  In  many  cases  the  solution  to  the  new  problem  differs 
from  the  body-alone  problem  only  in  interpretation. 

The  purpose  of  this  report  is  to  study  minimum  wave -drag  combina¬ 
tions  which  satisfy  a  few  of  the  many  possible  geometric  restraints  per¬ 
tinent  to  the  interests  of  airplane  designers.  An  attempt  has  been  made 
to  analyze  the  various  problems  in  a  unified  manner  so  that  extensions 
to  other  kinds  of  restraints  can  be  deduced. 
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LIST  OF  IMPORTANT  SYMBOLS 

aspect  ratio 

source  distribution  equivalent  to  wing  in  sense  defined  by 
equation  (3) 

multipole  distribution  of  order  n 
wave  drag 

portion  of  drag  due  to  all  the  nth  order  multipoles  for  n>0 
See  equation  (8). 

additional  drag  resulting  from  restraint  (See  eq.  (ll).) 
restraints  defined  in  equations  (19) 

distance  between  apexes  on  x  axis  of  forecone  and  aftercone 
enclosing  wing  (See  sketch  (c).) 

length  of  basic  body 

length  of  modification  to  basic  body 

Mach  number 

PQU02 

2 

average  body  radius 

fuselage  area  in  cross'  section  normal  to  the  free  stream 

normal  (to  free  stream)  projection  of  wing  area  in  section  cut 
by  plane  xx  =  x  +  £ y1  cos  6  (See  sketch  (b).) 

speed  of  free  stream 
volume 

Cartesian  coordinate  system  (See  sketch  (a).) 

source  distribution  representing  the  fuselage  modifications 
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(3  '  n/M2.-  1 

6  polar  coordinate  (See  sketch  (a).) 

pQ  free -stream  density 

a  See  equation  (17)- 

cp  velocity  potential 

BASIC  THEORY  AND  ASSUMPTIONS 
Basic  Theory 


Many  of  the  discussions  and  derivations  contained  in  the  following 
are  carried  out  on  the  assumption  that  the  reader  is  familiar  with  the 
concepts  presented  in  reference  k  which  should  be  considered  as  a  first 
part  to  this  report.  In  particular,  A  z 

an  acquaintance  with  the  solutions  .X 

to  the  wave  equation  referred  to  as  s' 

"multipoles"  is  assumed,  together 

with  Hayes T  invariance  principle  and  j  0 

the  consequent  multipole  distribu-  / / 

tions  equivalent  to  a  wing  in  the  /  \ _ 1 

sense  that  both  induce  the  same  A 

momentum  spectrum  at  infinity.  ^  £ - A 

The  entire  analysis  used  herein  / 

is  based  on  the  assumptions  and  ideali-  VJ' 

zations  necessary  to  develop  the  X 

linearized  equation  for  the  velocity 
potential,  cpj  in  supersonic  flow, 

namely  1 

32(Pyv  -  <fW  -  q>7.7.'  =  0  (!)  ✓ 


where  (32  =  M2  -  1  and  the  reference 
coordinate  system1  is  shown  in 
sketch  (a) .  The  analysis  is  further 
restricted  to  the  solution  of  prob¬ 
lems  involving  a  given  uncambered 


1It  should  be  stressed  that  the 
stream  direction  (wind  axes) . _ 


_ Sketch  (a) _ 

axis  is  parallel  to  the  free- 
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and  untwisted  wing  mounted  centrally  on  a  vertically  symmetrical  fuse¬ 
lage,  the  entire  configuration  being  at  zero  angle  of  attack. 


Additional  Assumptions 


We  now  make  the  two  additional  assumptions:  one,  the  value  of  (3A, 
where  A  is  the  wing  aspect  ratio,  is  small;  and  two,  the  value  of  (3R/Le > 
where  R  is  the  average  body  radius  and  2Le  is  the  distance  along  the 
x  axis  in  which  the  multipole  strengths  differ  signif icantly  from  zero, 
is  small.  " 

One  can  evaluate  the  significance  of  these  assumptions  by  studying 
their  implications  relative  to  the  source  and  multipole  distributions 
used  to  simulate  the  wing  and  body.  Suppose,  for  example,  a  group  of 
multipoles  are  placed  along  the  body  center  line,  their  strengths,  Cn(x), 
being  fixed  by  the  condition  that  a  circular  cylinder  in  the  vicinity  of 
the  body  is  a  stream  surface  when  the  velocity  field  induced  by  these 
multipoles  is  combined  with  the  velocity  field  induced  by  the  source  * 

sheets  representing  the  given  wing.  With  the  assumptions  of  small  (3A 
and  PR/Le  mentioned  above,  the  Cnfs,  for  n  greater  than  0,  can  be 
shown  (see,  e.g.,  ref.  4)  to  have  a  negligible  effect  on  the  wave  drag.  ^ 

Hence,  all  the  multipoles  (for  n  >  0 )  that  combine  with  the  wing  to  make 
a  circular  cylinder  a  stream  surface  and  any  additional  multipoles  (for 
n  >  0)  added  to  make  the  body  have  mild  distortions  from  such  a  surface 
are  negligible  in  evaluating  the  wave  drag.  Therefore,  under  the  assump¬ 
tions  mentioned  above,  out  of  all  the  singularities  distributed  along  the 
body  axis,  it  is  necessary,  in  studying  the  wave  drag,  to  consider  only 
the  sources . 

With  the  restrictions  to  small  values  of  PR/Le  and  mild  body  dis¬ 
tortion  (see  Ward,  ref.  6,  for  a  discussion  of  orders  of  magnitude), 
slender -body  theory  can  be  used  to  calculate  the  body  shape,  and  on  the 
basis  of  this  theory  one  can  show  (see  ref.  4,  Appendix  B)  that  Sf(x), 
the  body  cross-sectional  area  measured  normal  to  the  free  stream,  is 
completely  determined  by  the  axial  source  distribution  alone.  Hence,  if 
only  the  exposed  panels  of  the  wing  are  used  to  calculate  the  Cnfs,  C0 
is  negligible  and  the  entire  axial  source  distribution  <xQ(x)  is  related 
to  the  geometrical  properties  of  the  body  by  the  relation 


ao(x)  =  U, 


dSf 


=  UoSf'OO 


(2) 
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The  Wing  Equivalent  Source  Distribution  and  the 
Optimum  Cancellation  Sources 

Let  the  given  wing  lie  in  the  zx  =  0  plane.  According  to  Hayes1 
theorem  (ref.  7),  the  wing  equivalent  source  distribution  [a0(x)]v  is 
obtained  by  accumulating  on  the  xx  axis ,  at  a  distance  x  from  the 
origin,  all  the  wing  sources  intercepted  by  the  line  xi  =  x  +  f3yxcos  0, 
and  then,  for  a  fixed  x,  averaging  these  values  as  0  varies  from  0 
to  2jt .  Thus,  using  thin-airfoil  theory  to  relate  the  planar  source  sheet 
to  wing  geometry,  one  finds 


~  ao(x) 


S^-  (x ,  0 )  d0 


where  &w'(x?0)  =  d/dx[Sw(x,0) ]  and 
S w(x,0)  is  the  normal  (to  the  x  axis) 
projection  of  the  wing  cross-sectional 
area  intercepted  by  the  plane2 
X]_  =  x  +  3yxcos  0  as  shown  in 
sketch  (b) .  Without  the  addition  of 
further  restraints,  the  optimum  source 
distribution  along  the  xx  axis  is 
that  which  just  cancels  the  wing 
equivalent  source  distribution.  Fur¬ 
ther,  this  can  be  interpreted  directly 
in  terms  of  both  fuselage  and  wing 
geometry  by  means  of  equations  (2) 
and  (3)-  Thus,  with  no  further 
restraints,  the  best  fuselage  shaping, 
for  a  wing-body  combination  satisfying 
the  assumptions  discussed  above,  satis¬ 
fies  the  equation 


*(x)  =  -  -L 

V  2rt 


Sv'(x,e)a.e  (4) 


and  has  any  reasonably  smooth  cross- 
sectional  form.  Notice  that  the  total 


S(x,0 )=  Normal  projection 
of  wing  area  along  A  A 


x+fycos# 


Sketch  (h) 


2The  true  oblique  plane  is  given  by  the  equation 

=  x  +  Pyxcos  9  +  pzxsin  9 

but,  to  be  consistent  with  the  assumptions  basic  to  linearized  theory, 
the  variation  with  Zi  is  neglected. 


6 


NACA  TN  3667 


volume  taken  out  of  the  fuselage  is  exactly  equal  to  the  total  volume  of 
the  exposed  portion  of  the  wing.  Hence,  the  total  volume  of  the  modified 
combination  is  the  same  as  that  of  the  original  smooth  cylinder . 


The  Drag 


The  total  wave  drag  of  a  system  can  be  expressed  in  terms  of  its 
actual  or  equivalent  multipole  distributions  as 


(5) 


where  Dn  is  the  drag  caused  by  the  nth  order  multipoles  an(x)  and  is 
given  by  the  equation  ' 


% 

q. 


r~2  f  4x1  f  dx2an^n+1^(xi)an^n+1^'(x2)ln  |  xi  -  x2  | ,  n 
“”U°  J-L0'  -L0 ' 


0,1,2.  . 
(6) 


where  an^  (x)  represents  (b/bx)  an(x) .  Under  the  assumptions 

00 

given  above,  the  magnitude  of  E  Dn  is  small.  Let  us  designate  it  by 
De,  so  that,  in  general. 


D  —  2Dq  4*  Dg 


On  the  other  hand,  the  total  wave  drag  of  a  system  composed  of  the 
combination  of  a  wing  and  a  body  can  also  be  written  symbolically  as 


D  =  Dw  +  +  Db  (8) 

where 

Dw  =  drag  of  the  wing  alone 
Dfc,  =  drag  of  the  body  alone 
Dwb  =  interference  drag 
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The  various  components  of  wave  drag  defined  in  equations  (7)  and  (8) 
help  one  to  evaluate  more  readily  the  drag  reductions  that  can  be  realized 
from' appropriate  fuselage  formations.  Thus,  if  the  fuselage  shape  satis¬ 
fies  equation  (4),  the  total  wave  drag  of  the  combination  under  the  assump¬ 
tions  that  pA  and  PR/Le  are  small  can  be  written  either  as 


D  =  D£ 

or  as 

D  =  Dw  -  ®b 


t9) 

(10) 


If,  in  finding  the  fuselage  shape, 

(a)  the  multipoles  representing  a  wing  and  a  body  flying  separately 
are  assumed  to  represent  the  same  wing  and  body  when  combined  (i.e., 
the  shape  fields  can  be  superimposed) , 

(b)  the  multipoles  representing  the  fuselage  are  equal  in  magnitude 
but  opposite  in  sign  to  the  wing  equivalent  multipoles, 

then  equation  (10)  holds  without  the  assumption  of  small  |3A  and  gR/Le - 

In  subsequent  problems  we  will  discuss  the  effects  on  the  wave  drag 
and  fuselage  area  distribution  of  adding  certain  additional  restraints  to 
the  body  geometry.  The  addition  of  such  restraints  may  or  may  not  change 
the  relation  given  by  equation  (10),  but  they  must  always  add  a  term  to 
equation  (9)  so  that 


D  =  Dg  +  Dre 

Dre  >  0 

j 

WINGS  CENTRALLY  MOUNTED  ON  SLENDER  QUASI -CYLINDERS 


(ID 


This  section  is  devoted  to  the  solution  of  two  problems  involving  a 
given  uncambered  and  untwisted  wing  mounted  centrally  at  zero  angle  of 
attack  on  a  tube  that  is  cylindrical  forward  of  some  point  ahead  of  the 
wing.  The  problems  are,  in  both  cases,  to  find  the  area  distribution  of 
the  fuselage  behind  the  cylindrical  portion  that  will  minimize  the  wave 
drag  of  the  combination. 
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Shortening  the  Fuselage 


Remembering  the  assumptions  listed  at  the  beginning  of  this  section, 
let  us  consider  the  following  problem: 

(■*•)  Given  a  wing  and  a  slender  fuselage  having  the  same  normal 
cross-sectional  area  in  all  planes  ahead  of  the  plane  x=  -LQ' 

(see  sketch  (c)),  what  is  the  optimum  fuselage  area  distribution 
behind  the  plane  x  =  -Lq 1  if  the  fuselage  must  end  at  the  plane 


y 
/  \ 

✓  \ 


Of  course,  if  Z0  >  L0  (i.e.,  the 
body  modification  can  extend  over  the 
entire  range  enclosed  by  the  forecone 
and  after cone  enclosing  the  wing) , 
the  solution  is  already  given  by  equa¬ 
tion  (k) .  Hence,  in  the  following, 

Zq  ^  Lq • 

/  For  simplicity  of  notations, 
let  a0(x)  represent  the  sources 
along  the  fuselage  center  line  and 
.  aQ(x)  represent  the  wing  equivalent 
,N  source  distribution.  Then,  according 
to  equation  (6) 


kixUr,2 


dx2[a0'(xi)  + 


\  / 
\  / 


\  (*i)  ]  [aQ  1  (x2)+a0  ^(xg)  ]Zn|xi-x2 


Sketch  (c) 


where  from  the  conditions  stated  in  the  problem  and  the  geometric  inter¬ 
pretation  to  the  fuselage  sources  given  by  equation  (2),  a  (x)  is  zero 
for  values  of  x  outside  the  interval3  -L0!  <  x  <  ZG. _ _ 

^It  is  necessary,  for  equation  (6)  to  be  valid,  that  Oo(-L 0T)  and 
^O'^o)  zero.  This  implies  that  ogq(x)  must  be  continuous  and  if  the 
body  shape  is  given  by  equation  (2),  this,  in  turn,  means  that  the  stream- 
wise  gradient  of  body  cross-sectional  area  must  be  continuous.  It  was 
pointed  out  in  reference  k  that  a0(-L0T)  and  a0(L0)  will  both  be  zero  if 
the  wing  has  no  blunt  edges  along  which  the  normal  component  of  the  free- 
stream  Mach  number  is  unity  or  greater. 
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Consider  now  a  variation  of  equation  (12)  for  a  fixed  a0(x)  in  the 
interval  -L0 T  <  x  <  L0  and  a  free  variation  of  a0(x)  in  the  sub¬ 
interval  -L0*  <  x  <  ZQ.  There  results 


Integrate  once  by  parts  with  respect  to  xi  (since  the  variations  6a0(-L0f) 
and  5a0(Z0)  must  be  zero).  Then,  by  the  fundamental  lemma  of  the  calculus 
of  variations,  the  bracketed  term  must  be  zero  for  .  -LQ r  <  xx  <  lQ  and  one 
finds  the  condition 


a ’ (x2)dx2 

— - -  0; 

xi  -  x2 


<  Xi  <  l0 


(13) 


Equation  (13)  is  an  integral  equation  which  can  he  inverted  (by  methods 
such  as  those  outlined  in  ref.  8).  Inverting,  integrating  once  with 
respect  to  x,  and  expressing  Oo(x)  and  a0(x)  by  means  of  equations  (2) 
and  (3),  respectively,  one  finds 


Sf’(x)  =  - 


1_ 

2it 


2Jt 


(Ik) 


which  gives  the  optimum  fuselage  area  distribution  under  the  conditions 
and  assumptions  posed. 

The  wave  drag  of  the  combination  represented  by. equation  (l4)  can  be 
expressed  either  in  the  terms  defined  in  equation  (8)  or  (ll) .  Let  us 
first  consider  the  form  given  by  equation  (8) .  If  the  expression  for.  the 
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or,  alternately, 


Are  -1  r  °  r  ° 


l0  l0 


In 


xi  -  x2|  (l0  +  L0') 


[7 (xi  +  L0')(x2  -  Iq)  +7 (x2  +  L0')(xi  -  lQ)  ] 


-  dx!dx2  (15b) 


Constrained  Fuselage  Areas 


Another  class  of  problems  is  that 
in  which  the  magnitude  of  the  fuselage 
area  is  fixed  at  various  points.  Sup¬ 
pose,  for  example,  that  a  fuselage 
shaped  according  to  equation  (k)  had 
in  some  region  a  cross-sectional  area 
too  small  to  be  acceptable  for  some 
practical  purpose.  The  question  is, 
then,  what\is  'the  best  shape  for  given 
values  of  ininimum  fuselage  cross- 
section  area  at  given  planes  and  what 
is  the  penalty  in  wave  drag  caused  by 
the  added  constraints?  Before  consid- 
ering  the  general  case  of  an  arbitrary 
number  of  restraints ,  let  us  first  con¬ 
sider  the  simple  problem: 


(ii)  Given  a  wing,  what  (under  the  various  assumptions  given 
above)  is  the  area  distribution  of  the  adjoining  fuselage  which 
has  a  prescribed  area  at  three  given  stations  (the  initial,  the 
final,  and  an  intermediate  station  x  =  d x,  see  sketch  (d) )  and 
yields  a  minimum  wave  drag  for  the  combination? 


As  before,  let  a0(x)  represent  the  wing  equivalent  source  distribu- 
^^on*  Then  the  drag  caused  by  the  restraints  can  be  written 

Ln  T 


p  ♦  p 

~q~  =  2iru02  J  [ao(xi) +U0Sf'(xi)  ]dxi  J 

“Lq  -L, 


ao’(x2) +  U0Sf" (x2) 


Xi  -  x2 


dx2 

(16) 
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where  Sf(x)  is  the  unknown  fuselage  area  to  he  optimized.  For  simplic¬ 
ity,  replace  the  unknown  Sf(x)  by  cr(x)  where 


X  2rt 

o (x)  =  f  falii  d6  + 


-L, 


U, 


Sf(x)  =  ^  J  Sw(x,0)d0  +  Sf  (x) 


^("I'o)  -  S0 

ff(Lo)  =  S2 


}  (17) 


Let 


0O(X)  =  d(x)j  -Lp  <  x  <  di 

t^i(x)  =  cr(x) ;  dx  <  x  <  L0 
and  the  restraints  on  the  fuselage  area  give  the  relations 

=  J, 


(18) 


fdl 

23T 

p 

/  ff0'(x)dx  = 

Si  -  So  +  —  /  Sv(dx,0)de 

-L0 

o 


(19a) 


rL° 

r  2Jt  n 

p 

/  bx^xjdx  = 

s2  -  s  -  —  /  Sw(dx,0)d0 

d-1 

2n  J 

0 

=  Ji  -  J0  (19b) 

where  J0  and  Ja  are  constants  fixed  by  the-  given  constraints.  Notice 

Jx  =  (S2  -  S0)  (19C) 


so  the  constant  Jx  is  a  measure  of  the  difference  between  the  initial 
and  final  areas. 

Using  the  usual  variational  techniques,  we  can  write,  for  the  quan¬ 
tity  to  be  minimized, ■ 

di 
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Co'CxiJdx! 


dl  g0”  (x2)dx2  +  rL° 
Xi  -x2  J 


Oj."  (x2)dx2 
Xi  -  x2 


+  Aq  + 


°  ffi'CxJdx!  F  1  gQ"  +  f  °  Sll 


a  itt(x2)dx2  +  ^ 


xi-x2 


Xl  -X2 


By  taking  the  variation  and  satisfying  the  conditions  at  the  end  points, 
one  obtains  the  two  simultaneous  integral  equations 


ao"(x2)dx2  _  P  °  ffi" (x2)dx; 


-Lp  <  X!  <  d; 


Xl  -x2 


Xl  -x2 


f  ^ 

J  x-1  - 


(x2)dX; 


Xl  -  X2 


di" (x2)dx2 
Xi  -  x2 


di  <  xx  <  L< 


The  set  of  equations  (20)  is  identical  to  that  analyzed  by  Adams 
(ref.  9)>  page  I**-,  for  bodies  of  revolution  with  fixed  areas  at  the 
initial,  final,  and  an  intermediate  section.  In  the  interest  of  subse¬ 
quent  generalization,  however,  we  will  consider  its  solution  in  the 
following  way:  First  write  the  equations  (20 )  in  the  equivalent  form 


r  0  g” (x2 

J  Xi  — 


(x2)dX; 


Xl  -x2 


'  O’  >  -Ln  <  Xl  <  d; 


-  — ;  di<xi<L( 
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One  can  show  that 


A  +  Bx2  ±  L0  -  dxx2 

a  (x2)  =  7"- -  Cicosh  -2 - — j 

vL02  -  x22  Lo  x2  “  dx  I 


is  the  solution  to  the  integral  equation  (where  A ,  B,  and  Cx  are  con¬ 
stants)  since 


L0  B  -  Cxcos  j—  )  ;  -L0  <  xx  <  dx 

a"  (x2)dx2  J  L  WJ 

Xl-  x2  "A 


-at  B+Cxcos-M— J  ;  dx  <  Xx  <  L( 


which  satisfies  equation  (21) .  The  constraints  can  now  be  satisfied  by 
means  of  the  equations 


a'(x)  =  J  a"  (xx)dxx  sb  A  cos-1^  -  bJ  L02  -  x2  + 
-L0 


Ci  (di  -  x)cosh 


-1  Lp  -  dqx 
Ln  x  -  di  I 


fL02  -  d.!2  cos-1L^ 


a(x)  -  SD  =  A  x  cos_^2  )  +  n/l02  - 

VLo/ 


x2  -  2  x  vL02  -  x2  +  L0  cos-1 (  ~ 


g)"- 


^[(d1-x)2cosh-1-g^'dlX| 
2  L  L0  x-  di 


r(L02-di2 )  (Lq2  -  x2) 


rL02  -  di2(di  -  2x)cos“1(f^J  + 
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Since  cr'(-L0)  =  cx(L0)  =  0 


Solving  for  Ci  and.  B,  we  find. 


Prom  equation  (25)  the  fuselage  cross-sectional  area  can  be  written 
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"the  additional  specification  is  made  that  the  initial  and  final 
areas  are  the  same,  the  solution  simplifies  considerably,  since,  for  such 
cases,  Ji  =  0  and  equations  (2 6)  and  (27)  reduce  to 


Sf(x)  =  S0  - 


-jQ  M*,a)da  +  — 


(di,0)d0 


(Lc^dxxWao'-x^tLo2-^2)  -  L/tx-dO^osh-iislllh 

L0  |x  -  dx 


D  Db  2jtL0 


*  *  1  (Lq  -  di  ) 


2  „  g,2  (Si  -  So)  Sx-  So  +  /  Sw(d1,e)d0  (29a) 


(V2  -  di2) 


Si  -  S0  +  tt“  /  Sw(d: 


,0)d0 


(29h) 


Often  the  exact  statement  of  the  restraint  is  that  S(x)  shall  not 
be  less  than  Si  at  x  =  di-  In  such  cases  care  must  be  used  in  applying 
equations  (26)  and  (27)  or  (28)  and  (29),  since  they  are  only  valid  when 
the  fuselage  cross-sectional  area  at  di  is  exactly  Si.  If  such  is  the 
case,  equations  (26)  and  (28)  give  the  optimum  body  shape  only  if  Jr,  >  0 
that  is,  only  if  0  -  ’ 


Si  >  S0  - 


Sw(di,0)d0 


Otherwise  the  optimum  variation  of  area  is  given  by  equation  (4) . 
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Next  let  us  generalize  the  analysis  leading  to  equations  (26)  and  (27) 
by  considering  the  following  problem: 

(iii)  Given  a  wing,  what  is  the  area  distribution  of  the  adjoin¬ 
ing  fuselage  which  has  prescribed  areas  at  n  +  1  stations  (includ¬ 
ing  the  initial  and  final  ones  fixed  by  the  Mach  forecone  and 
aftercone  enveloping  the  wing,  see  sketch  (e))  and  yields  a  mini¬ 
mum  wave  drag  for  the  combination? 


By  analogy  with  equation  (22),  the  Original 

integral  equation  for  a!t(x)  (where  a(x)  Modified  Cylinder 

is  defined  by  eq.  (17))  that  must  be  body  ^  — 

satisfied  for  a  minimum  wave  drag  can  be 

written  at  once  in  the  form  *  c  !  \  si' v 


r  °  cr"(x2)dx2 

/  v  =  AiJ  di  <  X!  <  di+1. 

r' 

\ 

\ 

1  1 

V\ 

1 

J  X!  -  X2 

I 

I 

J- 

i  —  0,1,  •  •  • ,  n 

where  d0  =  -L0  and  dn+1  =  L0.  The 
quantity 


\5*T  I  dn/ 


\ 

\/  »y  Are  os  S0, 

SnH>and 
T  wing  ail 
™  fixed. 

Sketch  (e) 


0"(x2) 


A  +  Bx2 

fr  2  1/2 


Lo  -  X2 


n 

-I01 


_  Ln  -  Xpd-? 

cosh”1  — 2. - 

L0|x2  -dj_ 


is  a  solution  to  equation  (30)  since  it  yields 
L  ^  u  t  _ 


a" (x2)dx2 
xi  -x2 


=  -3t  B  +^CjCQS~1^^^  -  ^ 


Cicos’1^^!)  j 

vjo J  _ 


dv-i  <  Xj  dv ,  v  —  1,2,  .  .  n  +  1  (32) 


in  which 


C-r  cos”1 


MCA  TF  3667 


Further,  it  is  apparent  from  equations  (2k)  and  (25)  that,  with  the  con¬ 
ditions  o'(-L0)  = a '(L0)  =  0 


tf'(x)  =  -B^Lq2-  x2  +  V  CiCdi  -  x)cosh~1— °  -  Xdl 

La  L0  X  dj_ 


Sf(x)  =  So  -  i  /  S^x,©) d©  -  |  x  L02  -  x2  +  L^cos"1!^* 


2„  u-i  Lo  "  xdi 


x)  cosh' 


Loix-dii  r± 


-n/ Ln^  -  d-t“|  d,cos--L(  ^  IWl2-  x2 


The  wave  drag  due  to  the  restraints  can  he  obtained  by  using  equa¬ 
tions  (32)  and  (16).  Thus 


a 

2  ®  (°o  “  an+i)  +  \  ^  Ci  Jtffi  - 


C0cos 


-xf^l 


-an+icos-1^ 


(35a) 


or  in  terms  of  the  components  defined  in  equation  (8) 


D  _  £w 

q  ~  q 


^  +  B(S0 


n 

sn+i)  +  ^  Si  jtSj[  - 


Socos'1^ 


Sn+iCOs-i(  — i 
V  -Ln 


(35b) 


where  a±  -  ^(di).  Notice  a0  =  SQ  and  cxn+1  =  Sn+1,  so  when  SG  =  Sn+1, 


fire  _  n  V'  , 

a  2  /  Ci(ai  -  S0) 


(35c) 
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or 

;  =  T'T  +  ”i°l(Sl  -3o)  (35d) 

i=i 

Finally,  using  the  known  values  of  Sf(x)  at  dv,  v  =  0,  1,  .  • 
n  +  1,  one  obtains  the  n  +  1,  simultaneous  equations 


2jt 


which  determine  the  n  +  1  constants  B,  Ci,  C2,  .  .  •  9  Cn-  These,  in 
turn,  fix  the  shape,  through  equation  (3^),  and  the  wave  drag,4  through 
equations  (35) • 

WINGS  CENTRALLY  MOUNTED  ON  SLENDER  CLOSED  BODIES  OF  REVOLUTION 


In  the  preceding  section  the  interference  between  the  central  portion 
of  the  airplane  and  its  nose  and  tail  regions  was  neglected.  In  this  por¬ 
tion  we  will  consider  the  entire  fuselage,  assuming,  first,  it  is  a  slender 
closed  body  and,  second,  it  can  be  calculated  in  the  presence  of  the  wing, 
using  the  same  postulates  given  in  the  previous  section  under  "Basic  Theory 
and  Assumptions 

4 

It  has  come  to  the  author's  attention,  through  a  private  communica¬ 
tion,  that  W.  T.  Lord  has  obtained  a  solution  similar  to  the  above  and  is 
using  it  to  calculate  the  drag  of  bodies  of  revolution  having  their  areas 
specified  at  a  given  number  of  stations.  Such  a  method  has  the  advantage 
of  giving  the  lower  bound  to  the  drag  of  bodies  whose  areas  have  been 
measured  at  a  discrete  number  of  places  and,  further,  of  giving  a  value 
representative  of  all  area  variations  in  the  vicinity  of  the  calculated 
optimum. 
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Unlimited  Indentation  Length,  Fixed  Volume 


Let  us  first  consider  the  question: 

(iv)  Given  the  wing,  body  length,  and  total  volume  of  the  com¬ 
bination,  what  is  the  area  distribution  of  the  body  which  yields 
a  minimum  wave  drag  if  the  apexes  of  the  Mach  forecone  and' 
aftercone  enclosing  the  wing  lie  within  the  body  (see  sketch  (f)) 
and  the  specified  volume  is  large  enough  for  the  body  to  be  real? 

This  problem  can  be  solved  in  a  simple  manner  by  means  of 
the . following  lemma  discovered  by  R.  T.  Jones,  using  methods 
similar  to  those  introduced  in  reference  10. 

Designate  the  closed  body  which  has  a  minimum  wave  drag 
for  a  fixed  volume  and  length  as  a  Sears -Haack  body.  Then  the 
total  drag  of  a  Sears -Haack  body  and  any  other  wing  or  (slender) 
body  entirely  within  the  fore  and  after  Mach  cones  with  apexes 
at  the  tail  and  nose  of  the  Sears-Haack  body,  respectively,  is 
given  by  the  equation 


/  2Vi  \ 

D  "  Esh  V1  +  Vsi  )  *  Dl  (37) 

where 

dSH  wave  drag  of  Sears-Haack  body  alone 

Di  wave  drag  of  second  body  alone 

^sh  volume  of  Sears-Haack  body 

V'i  volume  of  second  body 

A  proof  of  this  lemma  can  be  obtained  by  placing  the  Sears-Haack 
source  distribution  and  the  wing  equivalent  multipole  distributions  (or 
the  second  body's  equivalent  multipole  distribution)  in  equation  (6). 
Since  only  the  sources  interfere,  the  drag  can  be  written  in  the  form 


where  the  interference  term  has  been  integrated  by  parts  and  -L0',  L0, 
and  -lQ,  lQ  form  bounds  of  the  arbitrary  and  Sears-Haack  source 
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and  since  one  can  easily  show 

rL° 

Vi  =  -  —  J  x1a0l(xi)dxi 

-Lo' 

equation  (37)  follows  Immediately. 

Returning  now  to  problem  (iv) ,  we  see  that  its  solution  follows  from 
equation  (37)  and  the  solution  is,  in  fact,  simply  a  Sears-Haack  body 
having,  at  the  appropriate  place  relative  to  the  wing -body  juncture,  the 
additional  area  variation  specified  by  equation  (k) .  This  follows,  since, 
if  Di  represents  the  combined  drag  of  the  wing  and  indentation,  then 
the  combined  volume  of  the  wing  and  indentation,  is  zero.  Hence,  the 
minimum  value  of  D,  for  a  given  volume,  is  obtained  when  Dgjj  and  Di  are 
independently  minimized.  But  Dgjj  is  already  a  minimum  on  a  volume  basis 
and  Di  is  a  minimum  for  a  given  wing.  Notice  the  location  of  the  wing 
along  the  body  is  immaterial,  provided  the  required  indentation  can  be 
accommodated  by  the  fuselage. 
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Limited  Indentation  Length  On  Sears -Haack  Body, 

Fixed  Volume 

Consider,  next,  the  more  difficult  problem 

(v)  Given  a  wing  and  Sears -Haack  body  of  length  2 ZQ  (long  enough 
to  contain  the  apexes  of  the  fore  and  after  Mach  cones  enclosing 
the  wing) ,  what  modification  of  this  fuselage  within  the  length 
^1  +  ^1  (and  within  that  length  only,  see  sketch  (g))  minimizes 

the  total  wave  drag  for  a  given  total  volume? 

In  order  to  answer  this  ques¬ 
tion,  it  is  necessary  to  consider 
separately  two  cases;  namely,  the 
one  in  which  Zi 1  >  L0r  and  Zx  >  L0 
(i.e.,  the  portion  of  the  body  free 
for  variation  contains  the  apexes  of 
the  wing's  Mach  cone  envelope,  as 
shown  in  sketch  (f))  and  the  other 
in  which  the  preceding  conditions 
are  not  satisfied. 

First  consider  the  combination 
for  which  Zx  *  >  L0 '  and  lx  >  L0. 

The  wave  drag  of  such  a  combination 
can  always  be  calculated  using  equa¬ 
tion  (37)  wherein  DgH  is  the  wave 
drag  of  the  basic  Sears -Haack  body 
fixed  by  the  stationary  nose  and  tail 
portions,  Dx  is  the  combined  wave 
drag  of  the  wing  and  the  (as  yet 
unknown)  body  indentation,  and  Vx 
is  the  net  difference  in  volume 
between  Vgjj,  the  volume  of  the  Sears - 
Haack  body ,  and  the  final  volume  of 
the  complete  configuration.  Since 
.  the  basic  Sears -Haack  body  is  fixed 

Sketch  (g)  and  the  total  volume  is  given,  the 

entire  term  Dgjj[l  +  (2Vi/Vgjj)  ]  is  fixed  and  the  solution  to  the  problem 
is  obviously  that  for  which  the  wing  equivalent  sources  and  the  source 
simulating  the  body  indentation  combine  to  form  a  Sears -Haack  distribu¬ 
tion  in  the  interval  -Zif<  x  <  Zx. 

Using  equations  (2)  and  (3)  to  relate  the  wing  and  body  source  varia¬ 
tions  to  their  respective  areas,  we  find  the  fuselage  cross-sectional  area 
can  be  written  for  -ZQ  <  x  <  -Zi* 
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for  -li'  <  x  <  1 1 

2Jt 

Sf(x)  =  (l02  -  x2)3/2  -  i  f  Sw(x,0)d9  + 

,  3  to  2*  J0 


- 12.^--  [(Z!>  +  x)(lx  -  x)]3/2  (4lb) 

3*(l  +  li') 

and  for  Zi  <  x  <  Z0 

■  Sf  (x)  =  (l02  -  x2)3/2  (to) 

3^Z0 


The  total  wave  drag  of  the  wing  and  the  fuselage,  as  given  by  equa¬ 
tions  (41)  is  then 

D  8VSH2  /  2VA  128V!2  '  p£ 

<1  rtl04  V  +  W  +  nih  +  h')4  +  ^ 

where  D€  is  defined  by  equations  (5);  (6),  and  (7). 

Since,  as  we  have  been  assuming,  PA  is  small,  De  is  negligible, 
and  a  comparison  between  equations  (37)  and  (42)  shows  that  the  drag  of 
the  combination  formed  by  mounting  two  wing  panels  on  a' Sears-Haack  body 
can  be  reduced  without  a  change  in  the  total  volume  and  with  a  modifica¬ 
tion  limited  to  the  interval  -Zi!  <  x  <  Zi  by  the  difference  between 

the  drag  of  the  two  panels  flying  alone  and  the  drag  of  a  Sears-Haack 
body  having  a  length  equal  to  Zi1  +  Zi,  and  a  volume  equal  to  the  volume 

of  the  two  panels.  So  long  as  the  points  x  =  -Zx*  and  x  =  Zi  do  not 

lie  off  the  basic  body,  and  so  long  as  the  required  indentation  can  be 
accommodated,  this  result  is  independent  of  the  wing 's* fore-and-aft 
position. 

If  the  body  modification  is  limited  so  that  either  Zx?<  L0f  or 
Zi  <  Lo  (see  sketch  (g))  or  both,  the  above  results / do  not  apply,  since, 
in  such  cases,  the  second  body  -  in  the  sense  defined  above  -  cannot  be 
varied  for  x  between  -Zx*  and  -LQT  or  Lc  and  Zi  or  both,  and  its 
drag  cannot,  therefore,  be  reduced  to  that  of  an  equivalent  Sears-Haack 
body.  The  best  modification  in  this  case  can  be  calculated  from  the 
results  presented  in  the  material  immediately  following. 
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Limited  Indentation  On  Arbitrary  Body  - 
Fixed  Volume 


Consider  the  question 

(vi)  Given  a  wing,  a  body  length,  and  the  area  distribution  of 
the  fore-and-aft  portions  of  a  body,  what  is  the  variation  of 
area  along  the  intermediate  portion  of  the  body  which  yields  a 
minimum  wave  drag  for  a  fixed  total  volume? 

Again,  as  in  equation  (17),  let  cr(x)  represent  the  sum  of  the  sources 
representing  the  basic  body  and  the  wing  equivalent  source  distribution. 


2% 

a(x)  =  Sf(x)  +  1  J  Sw(x, 


e)de 


(^3) 


It  is  now  convenient,  however,  to  let  c(x)  be  a  fixed  function  in  the 
entire  interval  -lQ'  <  x  <  lQ,  see  sketch  (h) ,  and  let  the  body  modifica- 

tions  y  which  are  to  he  optimized  in 
the  interval  i  <  x  <  2i  he  repre¬ 
sented  hy  ASp(x)  which  has  the  end 
conditions 


4Sf  *(-li)  =  ASf  ’(h)  =  0 
ASfi-h)  =  ASf(li)  =  0  j 


m 


The  change  in  volume  caused  by  the 
body  modification,  AV,  is  given  by 


AV 


<-1 

=  -  J  xASf  *  (x)dx  (45) 

-Zi 
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The  usual  variational  procedure  leads  directly  to  the  integral 
equation 


j'  ASf" (x2)dx2 
-lx  xl“x2 


,  a"  (x2)dx2 


J  7  r  xx  -  x2 

-  or\ 


+  A0  +  A1x1;  -lx  <  x  <  (46) 


where  Ac  and  Ai,  are  fixed  by  the  conditions  given  in  equations  (44) 
and  (45).  Equation  (46)  is  similar  in  form  to  equation  (13)  and  its 
inversion  can  be  obtained  by  use  of  methods  similar  to  those  for  invert¬ 
ing  the  latter  equation.  Thus,  the  solution  to  equation  (46)  becomes 
for  -2  1  <  x  <  Zi 


ASfr(x)  =  -a*(x)  + 


2  p 


lx  -  X' 


n  -aT(xi)dxi 

lQ  1  (xi  -  x)  \ /  Xi2  -  Z  x‘ 


r  — +  fa.rgIs-;^Io.2(T+AV)i  +  Ei 

Ii  (xi  -  x)-J  Xi2  -  Zi2  Ji4  *-  -*  ll 


where 


■>Zl  xin0'(xx)dxi  p°  xinc'(xi)dxi 


Xl  -  l 


Xi2  -  l x2 


and  V  is  the  total  volume  of  the  wing  and  unmodified  fuselage,  that  is 
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Equation  (4-7)  integrates  to  give 


.  „  If  Do  is  the  ^ag  of  the  original  wing-body  combination  and  DAQ 

is  the  drag  of  a  body  of  revolution  having  the  same  normal  area  distribu¬ 
tion  as  the  modification,  then 


D  _  Dff  8AV  r  -1 

*1  q  "  q  +  2^V  +  +  Jl2lo  ~  2I2  '  (52a) 


On  the  other  hand  Dre  can  be  written 
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Limited  Indentation  On  Arbitrary  Body  - 
Fixed  Diameter 


As  a  final  example  in  this  section,  consider  the  question 

(vii)  Given  a  wing,  a  body  length,  and  the  area  distribution  of 
the  fore-and-aft  portions  of  a  body,  what  is  the  intermediate 
variation  of  fuselage  area  that  has  a  given  area  at  some  inter¬ 
mediate  station  Xi  =  dx  and  yields  a  minimum  wave  drag  for 
the  combination? 

Using  the  same  definition  for  a(x)  as  is  given  in  equation  (43),  and 
again  designating  the  area  modification  as  ASf(x) ,  one  can  apply  the 
same  methods  used  to  develop  equations  (21)  and  (46)  and  write  the  inte¬ 
gral  equation  for  ASf(x)  in  the  form 


where  A0  and  Ar  are  constants  whose  values  depend  upon  the  restraints. 
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where 


Cl  =  - — -  Usf  (di)  +  a(da)  -  H(d!)  -  -Zl-~  Cil  (Ildl  +  I0Z!2)]  (60) 

(U2-di2)2  L  «U2  J 

and  In  and  H(x)  are  defined  in  equations  (48)  and  (51),  respectively. 
The  drag  can  be  expressed  either  as 

D  Dg  DAq 

-  =  —  -  —  +  nCiASf  (di)  (6la) 

where  D^  is  again  the  drag  of  the  original  unmodified  combination  and 
Das  is  "the  drag  of  the  modification  alone,  or  as 


51  +  B2(*i2-dig)  +  1_ 

1  4jt  4it 


(6lb) 


where  D  and  Ci  are  defined  by  equations  ($4)  and  (60) ,  respectively, 
and  B  is  given  by 


Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Moffett  Field,  Calif.,  Aug.  4,  1955 
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